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UDC  621.396.96'06 

ALLOWANCE  FOR  ANTENNA  ERROR  IN  DETERMINING  RADAR  TARGET  SCATTERING  MATRIX 

Kiev  IZVESTIYA  VYSSHIKH  UCHEBNYKH  ZAVEDENIY :  RADIOELEKTRONIKA  in  Russian 
Vol  26  No  3,  Mar  83  (manuscript  received  after  revision  10  Mar  82)  pp  74-75 

[Article  by  B.A.  Atayants  and  V.V.  Yezerskiy] 

[Text]  The  target  scattering  matrix  [1] 
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in  any  polarization  basis 
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where  and  T  ^  —  basis  angle  of  orientation  and  ellipticity, 
pb  pb 

respectively  —  is  done  with  antennas  having  two  orthogonal  polarizations 
which  coincide  with  the  unit  vectors  of  the  basis  (2).  However,  because 
of  the  finite  accuracy  with  which  antennas  are  fabricated  their  actual 
polarization  differs  from  the  required  figure.  It  is  important  to  investigate 
the  connection  between  discrepancies  in  measurements  and  discrepancies  in 
the  polarization  parameters  of  the  antennas,  and  to  find  a  method  to  allow 
for  errors. 

The  voltage  at  the  input  of  the  receiving  antenna  can  be  found  from  the 
formula  [2] 


=  Asb, 


(3) 


where  the  ^  designates  transposition,  and  the  column  matrices  A  and  B 
describe  the  receiving  and  transmitting  antennas.  The  elements  of  matrices 
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A  and  B  are  the  projections  of  the  radio  wav;^  formed  by  the  antennas  in  the 
radiating  mode  onto  the  unit  vectors  and  of  polarization  basis  (2). 

In  the  general  case  the  errors  in  the  antenna  orientation  and  ellipticity 
angles  differ  for  different  polarizations*  We  shall  designate  the  actual 
antenna  polarizations  T^) ,  and  T^) ,  respectively, 

for  the  receiving  and  transmitting  antennas.  Then 

,  At;,  =  t^  -  T  ^  n/2,  At^  =  , 

—  ■'l^b »  A't^B  ''^pb*  =  'ta  —  '*l’pb~'  At^  =  xg  +  x^^  .  (4) 


Knowing  the  results  of  measuring  the  voltage  at  the  output  of  the  receiving 
antenna  ^12’  ^22  three  different  antenna  polarization 


combinations,  we  can  obpain  from  (3) 
respect  to  the  unknown  S^2>  ^22 


the  following  system  of  equations  with 
(assuming  ^22^^21^  * 


^11  ==  +  (^1^2  +  ^2^l)  •^12  +  ^2^2'^22»  ^12  ”*  +  (^1^2  +  ^2^0  Ki  ^  ^2^2^22* 

^22  “  '^11  +  (^1^2  +  ^2^1 )  *^12  +  ^2^2‘^22- 


(5) 


The  coefficients  of  ^22  tkis  system  can  be  found  by 

using  the  unitary  matrix  obtained  in  [3]  to  calculate  the  transformations 
of  the  elements  of  column  matrices  A  and  B: 


a[  =  cos  A^l)^  cos  At^  +  /  sin  AiJ)^  sin  {2r^^  At^), 

—  /  cos  A-ip^  sin  At^,  Oj  =  —  sin  A'i))]]j  cos  (At^  —  4-  /  cos  Aij)^  sin  At^^» 

=  cos  A^l)^  cos  At^  +  /  sin  Ai]?'^  sin  (Ax^  —  ^^b^’ 


(6) 


b^,  b^,  bp  b”  are  described  by  analogy  with  (6),  substituting  the 
subscript  B  for  A. 

•  •  • 

By  solving  (5)  for  the  unknowns  ^22’  obtain  the  following 

in  general  form: 
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These  expressions  are  not  valid  when  the  determinant  of  system  of  equations 
(5)  is  zero,  which  can  occur  in  one  of  these  three  situations: 


^2^1  — 


*1*2  ■ 


*2^'!  =  0, 


“1*2  ■ 


•  02^1  =  0. 


(9) 


It  is  easy  to  demonstrate  that  one  of  these  conditions  can  be  satisfied 
in  accordance  with  one  of  these  practical  cases:  1) 

both  polarizations  of  the  receiving  antenna  are  circular  and  coincide; 

2)  t'  =t''  =+45°,  i.e.,  both  polarizations  of  the  transmitting  antenna  are 

circular  and  coincide;  3)T’^T’g=+^45°  or  i.e.,  V^2  ^21 

measured  with  the  receiving  and  transmitting  antennas  having  coinciding 
circular  polarizations. 

It  is  clear  that  these  situations  contradict  the  original  statement  of 
the  problem  of  measuring  the  scattering  matrix  by  antennas  with  different 
polarizations,  and  indicate  that  two  of  the  three  equations  in  system  (5) 
become  linearly  dependent.  Therefore,  it  is  impossible  for  conditions  (9) 
to  be  satisfied  in  practice. 

Expressions (7)  and  (8)  thus  make  it  possibly  to_^calculate  the  scattering 
matrix  elements  in  any  polarization  basis  (e^^,  e^)  by  the  results  of  three 

measurements  of  the  complex  voltages  by  antennas  with  polarizations  which 
are  different  from  those  required. 

Let  us  look  at  three  cases  which  are  of  practical  importance. 

1.  At’^=  A4j’g=Af’g,  AT’g=  “A^”b^  the  antenna 

errors  are  the  same  for  both  polarizations  and  differ  for  the  receiving 
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and  transmitting  antennas.  It  then  follows  from  (5)  that 


(10) 


where  *  designates  complex  conjugation. 

2.  Alp'  =  Alp'  ,  =Ai|j"  ,  At'.=At'  ,  At".=At''  ,  i.e.,  the  same  antenna  is 
used  for  receiving  and  transmitting.  Then 


(11) 


3.  =At'’  ,  At’  =  -At”  =At’.^=  -Ax’*  ,  i.e.,  the  conditions 

^A^A  B  B  A  A  B  B 

of  the  first  two  cases  are  satisfied,  which  yields 


(12) 


Consequently,  the  linear  relationship  (7)  holds  between  the  measurements 

...  ... 

of  V^2’  ^22  scattering  matrix  elements  S^2>  S22  the 

polarization  basis  (e^^,  62).  The  coefficients  which  enter  into  expressions 
(7)  are  in  general  functions  of  the  eight  complex  numbers  a'^,  a'25  ^"2 


b’^,  b'2»  ^"2  they  are  expressed  through  the  antenna  polarization 

parameter  errors  (6). 


Expressions  (8)  can  be  simplified  in  certain  cases  by  reducing  the  number  of 
unknown  variables  to  four  (10),  (11),  or  even  to  two  (12). 


Obviously,  these  findings  allow  the  a  priori  known  measurement  antenna 
parameters  t’^,  ip”^,  T"g  and  the  results  of  the 

three  measurements  of  ^22  calculate  the  elements  of  the 

scattering  matrix  in  any  polarization  basis. 

This  analysis^ is  valid  for  single-position  radar  location,  i.e.,  when  the 
equality  S.„=S„-  is  satisfied.  Violation  of  this  condition  makes  it 

necessary  to  make  four  measurements  of  ^21  ^22  four 

combinations  of  polarizations  of  the  receiving  and  transmitting  antennas 
and,  accoringly,  to  solve  a  system  of  four  linear  equations  instead  of  (5) . 
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The  final  expression  for  the  sought 
than  the  three  in  (7). 


will  contain  four  terms,  rather 
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COMMUNICATIONS 


UDC  621.391.8 

CROSS  SPECTRAL-CORRELATION  PROCESSING  OF  SIGNALS  IN  DIFFERENT  ORTHOGONAL 
BASES 

Kiev  IZVESTIYA  VYSSHIKH  UCHEBNYKH  ZAVEDENIY:  RADIOELEKTRONIKA  in  Russian 
.  Vol  26^No  3,  Mar  83  (manuscript  received  after  revision  28  Apr  82)  pp  3-7 

[Article  by  N.V.  Bebikh  and  A. I.  Denisov] 

[Text]  This  study  investigates  aspects  of  cross  spectral- 
correlation  transformations  in  different  orthogonal  bases. 

A  method  is  examined  for  computing  the  correlation  function 
by  the  method  of  spectral-correlation  transformation  and  of 
computing  the  spectra  and  correlation  functions  of  certain 
signals ;  the  structural  diagram  of  a  spectral-correlation 
converter  is  presented  as  well. 

Systems  of  non- sinusoidal  functions  have  recently  been  used  effectively  in 
modern  radio  electronics.  These  include  systems  of  Walsh,  Haar  and  other 
functions  which  are  convenient  for  digital  processing  by  computer,  as  well 
as  systems  of  Laguerre,  Hermite  and  Karunen-Loew  functions  for  investigating 
transient  attenuating  processes,  etc.  The  choice  of  a  system  of  basis 
functions  is  determined  by  the  mean  square  approximation  error,  the  amount  of 
computational  labor  Involved,  the  noise  tolerance  of  the  system  of  basis 
functions  and  the  ease  of  hardware  implementation.  In  order  to  provide  an 
operational  analysis  of  signals  under  study  in  the  most  suitable  orthogonal 
basis  it  is  helpful  to  have  the  capability  to  switch  quickly  between  different 
systems  of  basis  functions. 

The  connection  between  spectral  coefficients  in  two  different  systems  of 
functions  is  demonstrated  in  general  form  in  [1]. 

Let  us  represent  the  investigated  signal  by  a  Walsh  series  and  define  the 
Walsh-Fourier  cross-transition  matrix  [M^l:  [M^]=l/N[Sp]x[W]  ,  where 

]  _  matrix  of  partial  sums  of  trigonometric  functions;  [W]  —  transposed 

Walsh  matrix,  with  dimension  NxN,  N=2’^,  where  n  order  of  system  of  Walsh 
functions. 

The  elements  of  matrix  [S^]  are  determined  by  the  following  formulas:  for 
the  odd  matrix  rows 
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sin  kxdx, 


2  simjN 


2n{m—\)}N 


where  k=l,2,...,  N/2  —  number  of  odd  row  in  matrix;  for  even  matrix  rows 


2Jim/Ar 


I 


COS  IxdXy  ’ 


where  k=l,2,...,  N/2  —  number  of  odd  row  of  matrix;  £=1,2,...,  N/2  --  number 
of  even  row  of  matrix;  m=l,2,...,  N  —  number  of  matrix  column;  N=2’^  -- 
nijmber  of  digitization  intervals. 

The  Fourier  spectrum  of  the  signal  in  question  can  then  be  computed  by  the 
formula 


N 

^Fj  =  ^ 
r-1 

/=i 


where  a^.  —  spectral  Fourier  coefficient;  —  spectral  Walsh  coefficient; 

yi  —  Walsh-Fourier  cross-transition  coefficients;  or,  in  matrix  form,  [A.p]  = 

ij 

=  [M^]x[A^],  where  [A^]  —  column  vector  of  spectral  Fourier  coefficient; 
[A^]  —  column  vector  of  spectral  Walsh  coefficient;  [1^]  —  Walsh-Fourier 
cross-transition  matrix. 


The  number  of  operations  required  to  make  the  cross  transition  between 
spectral  Walsh  and  Fourier  coefficients  can  be  defined  by  the  recursive 
formula 


Lwf  ~  1  +  2  ^ 

1 


(2"“*  -  1), 
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where  n=l,2,.,,  —  order  of  system  of  Walsh  functions;  N=2^  — number  of 
spectral  coefficients. 

The  Walsh-Haar  cross-transition  matrix  and  Harr-Fourier  matrix  [M^^] 

can  be  calculated  analogously: 


[Mwh]  =••  UN  [S//]  X  [Mhf]  —  UN  [§/:•]  X  [Hf, 


where  [S  ]  — 
n 

of  normalized 


T 

matrix  of  partial  sums  of  Harr  functions;  [H] 
Haar  functions  matrix,  with  dimension  NxN. 


transpose 


The  number  of  operations  required  for  the  Walsh-Haar  cross  transition  can 
be  calculated  by  the  formula 


The  number  of  operations  for  the  Haar-Fourier  cross  transition  is  defined 
as  L^p=2(L^_^-l)+(2’^-l-l)(2’^  -1),  n=2,3...,  for  n=l, 

The  total  number  of  operations  required  to  compute  the  spectral  coefficients 
in  Walsh,  Fourier  and  Haar  bases  with  the  help  of  the  cross-transition 
matrices  and  will  be:  =L^+L2+L2"^Ii^,  where  Lj^=2’^ 

computation  of  partial  sums  of  signal  in  question;  L2=n2’^  —  calculation  of 

Walsh  coefficients  with  help  of  fast  Walsh  transform;  —  calculation 

of  Fourier  coefficients  with  help  of  Walsh-Fourier  cross  transition  matrix; 

L,=L„,  —  calculation  of  Haar  coefficients  with  help  of  Walsh-Haar 
4  WH 

transition  matrix.  For  N=16  L  =16+64+71+156=307 . 

^°^W 

When  a  system  of  Haar  functions  is  used  as  the  initial  system,  the  total 
number  of  calculations  will  be  L.  .  =LT4'Lr>+LQ+h,  ,  where  L-=2^  calculation 

totH  12  3  4  1 

of  partial  sums  of  signal  in  question;  L2=2(2^-l)  —  calculation  of  Haar 
spectrimi  with  help  of  fast  Haar  transform;  13=!^^  —  calculation  of  Walsh 
coefficients  with  help  of  matrix  >  ^4~^HF  —  calculation  of  Fourier 

coefficients  by  matrix.  For  N=16  1^^^  =16+30+71+329=446. 

If  the  spectral  coefficients  are  calculated  in  Walsh,  Fourier  or  Haar 
bases  without  using  cross  transition  matrices  and  with  the  help  of  fast 
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transforms,  the  total  number  of  operations  will  be  L^^^=L^+L2+L2+L^.  where 

L^=2’^  —  calculation  of  partial  sums  of  signal  in  question;  L2^n2^ 

calculation  of  Fourier  coefficients  with  help  of  FFT  algorithm  (number  of 
complex  multiplications  followed  by  addition  and  subtraction) ;  L^— n2 

calculation  of  Walsh  coefficients  with  help  of  fast  Walsh  transform  (number 
of  additions  and  subtractions);  L^=2(2’^-l)  —  calculation  of  Haar 

coefficients  with  help  of  fast  Haar  transform  (number  of  additions  and 
subtractions) . 


The  use  of  cross  transition  matrices  to  calculate  spectra  in  different 
bases  reduces  the  amount  of  computation,  with  the  greatest  gain  being 
achieved  when  a  system  of  Walsh  functions  is  employed  as  the  initial  basis 
system.  When  implementing  the  device  in  hardware,  in  order  to  determine 
the  spectral  coefficients  in  different  orthogonal  bases  it  is  best  (in  terms 
of  hardware  costs)  to  determine  the  spectral  coefficients  in  one  basis 
(e.g.,  a  Walsh  function  basis),  and  then  to  switch  to  other  function  bases 
with  the  help  of  cross-transition  matrices,  which  will  consist  of  weighted 
resistance  units  and  adders  in  terms  of  hardware. 


Figure  1  shows  the  structural  diagram  of  the  spectrum  converter.  The 
spectral  coefficients  from  the  output  of  Walsh  spectrum  analyzer  1  are 
input  to  the  inputs  of  the  weight  coefficient  units  —  Fourier  3,  Haar  4 
and  Laguerre  5,  the  outputs  of  which  are  connected  to  switch  7.  A  signal 
from  control  system  8  causes  the  switch  to  connect  the  proper  weight 
coefficient  unit  to  adder  unit  9,  which  outputs  the  spectral  coefficients 
in  the  proper  orthogonal  basis  in  parallel.  Converter  10  is  used  to  obtain 
the  serial  spectrum. 

This  device  permits  spectral-correlation  conversion  of  signals.  In  order 
to  do  this,  synthesized  signal  distribution  unit  2  and  multiplier  unit  6 
must  also  be  used.  The  correlation  function  K(t)  in  this  case  is  computed 
as  follows : 


K{-t) 


N 


m=0 


j+m, 


where  m=0,l,2,...,  N-1;  j=l,2,3, . . . ,N;  —  value  of  synthesized  signal  in 

digitization  interval;  T=m/N  —  shift  parameter. 

Since  Walsh  functions  are  piecewise  constant  over  the  digitization  Interval 
1/N,  it  is  sufficient  to  determine  the  values  of  the  correlation  function 
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at  m/N  points,  assuming  the  changes  between  those  points  to  be  linear  [2]. 
Distribution  unit  2  connects  to  the  inputs  of  multipliers  6  the  required 
values  of  the  synthesized  signal,  which  are  known  in  each  digitization 
interval,  as  a  function  of  the  value  of  the  shift  parameter  T.  The  order 
of  the  switching  is  assigned  by  converter  control  system  8,  With  this  method 
of  obtaining  the  correlation  functions  there  is  no  need  to  provide  a  time 
offset  T  between  the  multiplied  signals  (like  in  ordinary  correlators,  where, 
e*g,  frequency-dependent  delay  lines  are  used  for  this).  In  the  present  case 
the  time  offset  is  assigned  by  simply  changing  the  combination  of  pairs  of 
factors  from  the  values  of  the  synthesized  signal  known  in  each  digitization 
interval. 

When  implementing  spectral-correlation  transformations  by  computer,  it  is 
best  to  represent  the  expression  for  the  correlation  function  as  follows: 


r 

W  =  ^  j*  = 

0 


dj  j*  Wal  (/,  0  Wal  (i,  t  +  t)  df, 

0 


where  a.  —  spectral  coefficients  of  expansion  of  signal  x(t)  into  Walsh 
1 

series;  a.  —  spectral  coefficients  of  expanding  the  signal  y(t)  into  a 
Walsh  series. 

The  internal  integral  represents  the  cross-correlation  function  of  the 
actual  Walsh  functions  and  is  independent  of  the  type  of  signals 

being  investigated;  therefore,  it  can  be  calculated  in  advance  for 
different  values  of  x. 

Then  the  expression  for  the  correlation  function  can  be  obtained  in  this 
form: 

/=i  1=1 


If  we  have  the  pre-calculated  matrices  [K^(t)]  for  various  values  of  T 

and  the  column  vectors  of  the  spectral  coefficients  of  the  signals  in 
question  [A^] ,  calculate  the  correlation  function  by  the 

matrix  method  in  accordance  with  this  algorithm: 

1.  Calculation  of  column  vector  of  inner  sum:  [cr^(T) ]=[K^(t^) ]x[A^] . 

2.  Formation  of  auxiliary  matrix  [cf(T)] 
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[a  (t)1  = 


ai(Ti)  ojTi)  ...cr^C-ti)  ■' 

*^1(^2)  ^2(^2)  •••^^(^2) 

L  ("^m)  0^2  ("'m)  ■•■^J  i'^m)  J 


where  j=l,2,...,N  —  number  of  matrix  row;  m=l,2,.,.,N  —  number  of  matrix 
CO lumn . 


Fig.  1 


3.  Determination  of  correlation  function  matrix 

[K  (t)]- [a  (T)]X  [>!;]• 

In  the  present  case  the  matrices  [K^(t)]  are  calculated  directly  during 

computation  of  the  correlation  function;  the  intermediate  results  of 
computing  [K^(t) ]  for  different  values  of  T  are  not  stored  in  memory,  but 

are  used  immediately  to  calculate  the  current  values  of  the  inner  sums  of 
[cr.(T)],  and  each  subsequent  matrix  [K^(t)]  is  written  in  the  location  of 

the  preceding  one,  which  reduces  the  computer  memory  size. 
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A  program  for  calculating  the  correlation  function  by  the  spectral- 
correlation  transformation  method  has  been  written  in  ”Assembler-BASS" 
language  and  has  been  implemented  on  an  ”Elektronika  NTs03-T”  microcomputer 
for  N=16.  A  standard  matrix  multiplication  program  and  matrix  row 
calculation  subprogram  were  employed. 

The  results  of  calculating  the  spectral-correlation  characteristics  of 
some  functions  in  Walsh,  Haar  and  Fourier  bases  are  shown  in  the  form  of 
time  diagrams  in  Fig.  2,  where  S(t) —  investigated  signal;  a^,  a^,  a^ 

Walsh,  Haar  and  Fourier  spectra,  respectively;  S’ (t)  —  synthesized  signal; 
K(x)  —  correlation  function. 

Analysis  of  the  findings  makes  it  possible  to  assess  the  effectiveness  of 
using  a  particular  function  basis  for  signal  processing  (by  minimum  mean 
square  error,  number  of  non-zero  expansion  coefficients,  correlation 
coefficient,  etc.,)i.e.,  it  permits  adaptation  to  the  signals  being 
investigated.  This  increases  the  effectiveness  of  quick  analysis  of 
complex  signals  and  of  synthesizing  optimal  correcting  and  control  effects 
in  automatic  information  processing  and  control  systems. 
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[Text]  This  study  solves  the  problem  of  computing  the 
power  spectrum  of  a  pulsed  random  process  with  determinate 
timing  Intervals  in  finite  closed  form.  A  computer  program 
is  described  for  calculating  the  power  spectra  of  signals 
controlled  by  a  regular  Markov  chain.  A  decomposition 
method  is  proposed  which  can  be  used  to  compute  the  spectra 
of  digital  signals  formed  from  packets  of  varying  length. 

Two-level  digital  signals  with  an  assigned  timing  frequency  are  employed 
widely  in  modern  communications,  hydro- acoustic  and  computer  equipment  [1]. 

When  deciding  which  of  the  existing  digital  signals  is  most  suitable  for 
effective  transmission  through  a  given  propagation  medium,  it  is  important 
to  be  able  to  find  the  power  spectrum  of  the  signal.  Knowledge  of  the 
power  spectrum  makes  it  possible  to  determine  the  amount  of  influence 
between  bits  within  symbols. 

The  limited  bandwidth  of  a  real  data  transmission  channel  causes  the 
occurrence  of  high  frequency  and  low  frequency  signal  distortions.  These 
distortions  can  be  reduced  by  selecting  signals  which  have  the  desired 
spectrum  shape  [2]. 

Two-level  digital  signals,  which  are  mappings  of  fixed-length  binary 
sequences  (block  codes) ,  have  attracted  a  great  deal  of  attention 
recently  in  digital  data  transmission  tasks.  Each  symbol  in  the  input 
binary  sequence  is  transmitted  by  a  pulse  with  duration  T,  or  by  a 
group  of  pulses  (a  "packet").  In  this  case  the  sequence  of  states 

of  the  digital  signal  is  modeled  well  by  a  process  controlled  by  a 
finite  ergodic  Markov  chain  [3].  Signals  formed  from  "packets"  of 
different  lengths  are  also  important  for  practical  purposes.  The  study 
[3]  solves  the  problem  of  computing  the  power  spectra  of  signals  formed 
from  "packets"  of  the  same  length  which  are  used,  e.g. ,  in  digital  magnetic 
recording  technology. 
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The  present  study  proposes  a  method  for  solving  the  problem  of  calculating 
the  power  spectrum  when  the  "packets"  forming  the  signal  are  of  different 
lengths. 


Formulas  for  calculating  power  spectrum.  The  pulsed  random  process  in 
question  consists  of  "packets"  with  duration  rT,  and  the  forms  of  the 
"packets"  are  defined  ^2y’***’  process  is  controlled 

by  a  finite  ergodic  Markov  chain  with  period  d  and  transition  probability 
matrix  P.  We  shall  designate  A=(A^,A2 .  A^)  the  probability  vector 

characterizing^the  stationary  distribution  of  the  Markov  chain  in  question; 


Aq  =  lim  P^’^,  ps=pds_A^^  s=l,2,...;  P^=I  —  identity  matrix; 
n-H» 

A=(A  6  ).  .  1  ;  <5..  —  Kronecker’s  symbol;  m^^’^^-P^P^e^Ae  ;  y=0>1j' 

1  ij'i,j=l,...,q’  13  T,  TJV  V 


vO_i 


(y)  Y  ^  (0)  .  ?  ?  %  -iwt 

with  duration  T  in  "packet";  q  —  number  of  "packets". 


dt;  r  —  number  of  pulses 


The  formula  derived  in  [3]  for  the  averaged  power  spectrum  of  a  random 
pulsed  signal  with  "packets"  of  the  same  length  controlled  by  an 
arbitrary  finite  ergodic  Markov  chain  can  be  rewritten  in  the  form 


P((0): 


4|ai^ 

rT 


d-l 


SIS' 

ji,V=l  Y=0 
00 


[Re  (I  [/  —  (P“  —  Ao)  + 


2  ^  —  23x/i/dr7)  /lo j  ^ 


h=—co 


sin  (dr/2  —  yr  —  p  +  v)  cor  .  T  , , - ,A 


(1) 


where  [I-(P^-AQ)e^'^’^^^]  ^  —  the  inverse  matrix  of  I-(P^-AQ)e^^’‘^'''. 

Formula  (1)  provides  the  solution  to  the  problem  of  calculating  the  power 
spectrum  in  finite  closed  form.  It  is  useful  in  that  it  allows  standard 
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programs  for  calculating  the  inverse  matrix  to  be  employed  in  computing 
the  power  spectrum. 

The  next  section  examines  a  numerical  algorithm  for  calculating  the  spectrum 
for  the  case  of  a  random  pulsed  signal  controlled  by  a  regular  Markov  chain, 
which  is  the  most  important  case  in  practical  terms.  In  this  case  the 
formula  appears  as 


/>(C0): 


4|a  M  1 


rT  2 


m^v  cos  ([X  —  v)  0)7  + 


ji.v—l 


+  y  y  mJi'Uos  (sr  +  p  —  v)  0)7 + 

s=l 


r 

+  COS  (|X  —  v)  ©r 

H.v=l 


2  ^  rT 


00 

h—~co 


2nh 


rT 


(2) 


<7 

We  recall  that  the  discrete  component  in  (2)  is  lacking  if  0,  VM- 

/=! 

In  particular,  this  is  the  case  when  the  matrix  P  is  doubly  stochastic 

Q 

and  2e/n  =  0,  VH- 
/=! 

Algorithm  for  calculating  power  spectra  of  signals  controlled  by  regular 
Harkov  chain.  Figure  1  shows  the  flowchart  of  the  numerical  algorithm  for 
computing  the  spectrum  in  accordance  with  formula  (2).  The  following 
notation  is  employed: 

r  ^ 

Q,  ==  -i  V  cos  (p  —  v)  (oT;  ^  ^  cos  (sr  +  p  —  v)  ©T ; 

H.V=1  s=l 


r 

Qs  =  ^  /nw  cos  (p  —  V)  ©r  f- 


00 
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The  following  is  essential  in  the  use  of  this  algorithm.  By  using 
corollary  (4.1.5.)  from  [4],  we  can  demonstrate  that  the  maximum  number 
of  required  iterations  needed  to  ensure  the  required  accuracy  A  of  the 

spectrum  computations  is  determined  by  the  formula 


SM  =  A^[lnA/ln(l  —28)  +  1], 


(3) 


where  N  —  power  of  matrix  of  transitional  probabilities  P  in  which  all 
elements  are  positive;  e  —  smallest  element  of  matrix  (actually,  the 
number  of  iterations  is  always  much  smaller) . 


Key: 

1  —  start 

2  —  input  [quantity] 

4  —  formation  of  system 

13  - store  [quantity] 

14  —  output  [quantity] 


Fig.  1 

15  —  if  S  =1,  =0;  if  S^^l  KS^  <S^ 

M  ’  ts  M  ts—  M 

16  —  input 

^  ts 

22  —  store  [quantity] 

23  —  output  [quantity] 

(continued) 
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Key  to  Fig.  1  (continued) 

24  —  Is  it  the  case  that  '5j^-S^3=0;  29  —  +Q3  “/“t 

if  not,  increase  S  ^  31  —  output  [quantity] 

ts  32  __ 

25  —  input 

26  —  output  [quantity] 


After  S,M  is  determined,  actions  are  taken  to  find  the  final  probability 

vector  A,  i.e.,  the  probabilistic  solution  of  the  algebraic  system  AP=A 
is  found  (e.g.,  by  the  Gauss  method  with  selection  of  i^he  main  element). 

We  note  that  if  P  is  a  doubly  stochastic  matrix,  then  A=(l,l,...,  l)/q. 

S=1  is  then  assigned  and  the  second  term  Q2  of  formula  (2)  is  computed. 

If  s  =S=1,  Q„  is  stored  and  then  and  Q„  are  calculated,  followed  by 
M2  X  j 

the  spectrum  P(to/(i)^)  for  all  given  w/o)^. 

~S+1 

If  S  ^1.  an  exhaustive  search  1<S  <S  is  made.  P  =P  P  is  first  computed, 
’  ts—  m 

and  then  the  procedure  described  above  is  repeated  cyclically. 

If  the  signal  in  question  has  a  large  transition  probability  matrix,  'ts 
may  be  large.  In  most  practical  cases  the  quantity  can  be  bounded. 

Then  the  algorithm  in  question  can  be  used  to  calculate  the  maximum  error 

6  for  a  given  S  in  accordance  with  the  formula 
II  ®  ts 


where  Q?  —  value  of  Q„  in  previous  iteration.  If  the  error  is  not 

excessive,  S  is  increased  no  further, 
ts 

In  the  below  example  of  calculating  the  spectrum  by  the  algorithm  described 

above,  the  null  members  of  the  matrix  P  include  those  elements  whose  values 

do  not  exceed  10~^.  For  S^^=l  the  calculation  took  1.33  minutes,  and  for 

S  =60  —  12.06  minutes, 
ts 

Decomposition  method  in  computing  power  spectra  of  signals  formed  from 
"packets"  of  different  lengths.  The  above  referred  to  the  case  in  which  the 
"packets"  from  which  the  signal  was  formed  were  of  the  same  length  rT.  A 
whole  group  of  signals  which  are  promising  for  use  in  digital  magnetic 
recording  and  which  are  obtained  by  means  of  group  recoding  of  the  input 
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binary  information  in  which  a  particular  number  of  zeros  or  ones  may  not 
occur  in  a  row  [5,6]  do  not  have  this  property.  These  signals  include  the 
(0,l)~code  or  Gabor  code  [5],  the  (4, 9) -code  or  Franaszek  code,  the  M^FM 
[6]  and  ChMP-1  signals  (Fig.  2,  a-d) . 


Ft 


^3/^32\ 


^3/^52^53 


(4,^)— codie. 


/ 


0 


0 


F 


0 


0 


0 


0 


3Tr 


b) 

Fig.  2 


$2 

3s/3s2 


3^/ 3^2 


0 


/ 


2rr 


c) 


ChMP-1 


0 

O/ 

/ 

0 

^2^22 

_J 

/ 

71 

If 

d) 

The  power  spectra  of  these  signals  are  successfully  calculated  by  the 
decomposition  method,  which  allows  these  signals  to  be  put  into 
correspondence  with  a  certain  Markov  chain  which  controls  the  decomposed 
random  process.  The  possibility  of  this  decomposition  follows  from  the 
fact  that  the  causal  connections  which  determine  the  structure  of  the 
process  for  the  class  of  pulsed  random  processes  with  determinate  timing 
intervals  in  question  correspond  precisely  to  the  probabilistic  connection 
between  its  elements. 

The  "packets”  from  which  the  signals  in  question  are  formed  are  thus 
divided  into  a  definite  number  of  "subpackets"  of  equal  length,  the 
probability  of  occurrence  of  which  as  part  of  its  own  "packet"  is  unity. 
The  dimensions  of  the  transition  probability  matrices,  of  course,  are 
greater. 

The  application  of  the  decomposition  method  is  illustrated  by  the  example 
of  computing  the  averaged  power  spectrum  of  the  signal  whose  structure  is 
explained  in  Fig.  2b, 
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Example  of  calculating  power  spectrum.  The  (4,9)-code  (Fig.  2b)  is  a 
promising  one  for  use  in  high  density  digital  magnetic  recording.  This 
signal  has  a  rather  complex  structure.  By  using  the  decomposition  method 
a  (4,9)-code  can  be  represented  by  24  alternating  "packets".  The  "packet" 
length  is  T  ,  r=6  (T  =6T) ,  q=24,  and  the  transition  probability  matrix 


The  matrices  e.-e, ,  P- ,  P»  are  shown  in  Fig.  3. 
1  o  1  Z 
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Fig.  3 


It  is  easy  to  note  that  a  transition,  e.g.,  from  "subpacket"  to 

"subpacket"  S-„,  has  a  probability  of  occurrence  of  unity,  since  both  of 
them  are  part  of  the  "packet"  S^,  while  the  probability  of  a  transition 
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from  ”subpacket”  (which  is  part  of  "packet"  S^)  to  "subpacket" 

(which  is  part  of  "packet"  S^)  is  1/8,  We  are  easily  convinced  that  the 

pulsed  random  process  which  is  obtained  is  controlled  by  a  regular  Markov 
chain.  In  order  to  use  formula  (1)  on  a  computer,  the  values  of  A  and  03/0)^ 

must  be  added  to  the  data  above.  We  shall  assign  A=10“^,  and  co/o)  =0-6  in 
steps  of  0.1. 

The  results  of  calculating  the  averaged  power  spectra  are  shown  as  plots 
in  Fig.  4,  where  the  values  of  P(a)),  multiplied  by  the  coefficient  7r^/2T^ 

are  laid  off  along  the  ordinate,  and  the  values  of  the  normalized  frequency 
are  laid  off  along  the  abscissa. 


Fig.  4 


Discussion  of  results.  Conclusions.  It  is  apparent  from  the  plots  in 
Fig.  4  that  the  maximum  spectral  energy  of  the  (0,l)-code  (curve  a)  lies 
near  the  normalized  frequency  a)/(jo^=l,  i.e.,  in  the  region  of  the  center 

frequencies  of  the  passband  of  a  magnetic  record-play  back  channel.  The 
constant  component  in  the  spectrum  of  the  signal  is  small.  These 
characteristics  of  the  spectrum  allow  the  (0,1) -code  to  be  recommended  for 
use  when  the  recording  density  does  not  exceed  256  bits/mm. 
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The  (4,9)-code  (Fig.  4b)  has  a  very  useful  spectral  shape  from  the  viewpoint 
of  increasing  recording  density.  The  maximum  spectral  energy  of  the  (4,9)~ 
code  lies  in  the  vicinities  of  the  normalized  frequency  (jo/a3^=0. 25.  However, 

the  spectrum  of  the  (4, 9) -code  contains  a  constant  component  twice  as  large 
as  that  in  the  previous  case,  which  results  in  substantial  distortions  to 
the  signal  as  it  passes  through  the  record-play  back  channel. 

2 

The  M  FM  signal  (Fig.  4c)  differs  from  the  others  in  Fig.  4  in  that  there  is 

no  constant  component  at  all.  The  maximum  value  of  the  spectrum  of  the 

M^FM  signal  occurs  at  the  normalized  frequency  a)/a}^=0.3.  These  characteristics 

of  the  M^FM  signal  spectrum,  which  are  helpful  for  increasing  recording 
density,  have  resulted  in  its  increasing  use,  especially  for  use  with 
floppy  disks. 

2 

It  should  be  noted  that  the  spectrum  of  an  M  FM  signal  as  calculated  by  the 
author^ s  method  and  the  spectrum  of  the  same  signal  cited  in  [7]  are 
exactly  the  same. 

Curve  d  in  Fig.  4  reflects  the  spectrum  of  a  ChMP-1  signal.  The  substantial 
level  of  the  constant  and  extreme  low-frequency  components  of  the  ChMP-1 
signal,  which  are  not  passed  by  the  magnetic  record-play  back  channel, 
cause  significant  distortions.  The  smallest  time  intervals  of  a  ChMP-1 
signal  which  corresponds  to  data  zeros  are  practically  not  reproducible. 

This  makes  it  much  harder  to  use  the  ChMP-1  and  similar  signals  in  high 
density  digital  magnetic  recording. 

BIBLIOGRAPHY 

1.  Bylyanski,  P. ,  Ingrem,  D.  "Tsifrovyye  sistemy  peredachi"  [Digital 
Transmission  Systems].  Moscow,  Izdatel’stvo  ”Svyaz’",  1980. 

2.  Polyak,  L.M.  "Veroyatnostnyy  analiz  mnogourovnevykh  signalov" 
[Probabilistic  Analysis  of  Multilevel  Signals].  Moscow,  Izdatel’stvo 
”Svyaz’",  1980. 

3.  Pletneva,  T.G.,  Rozorinov,  G.N. ,  Eydel’man,  S.D.  "Averaged  Power 
Spectrum  of  Pulsed  Random  Processes  with  Determinate  Timing  Intervals". 

IZV.  VUZOV.  RADIOELEKTRONIKA,  1981,  Vol  24  No  8  pp  31-38. 

4.  Kemeni,  DZh. ,  Snell,  DZh.  "Konechnyye  tsepi  Markova"  [Finite  Markov 
Chains].  Moscow,  Izdatel’stvo  "Nauka",  1970. 

5.  Gabor,  A.  Adaptive  coding  for  self-clocking  recording.  IEEE  Trans. , 

1967,  V.  EC-16,  N  6,  pp  866-868. 


21 


6.  Mackintosh,  N.D.  The  choice  of  a  recording  code.  The  Radio  and 
Electronic  Engineer,  1980,  v.  50,  N  4,  pp  177-193. 

7.  Lindholm,  D.A.  Power  spectra  of  channel  codes  for  digital  magnetic 
recording.  IEEE  Trans.,  1978,  v.  MAG— 14,  pp  321-323. 

COPYRIGHT:  "Izvestiya  vysshikh  uchebnykh  zavedenly  MV  i  SSO  SSSR. 

Radioelektronika",  1983 

6900 

CSO:  8144/1267 


22 


UDC  621.391.14:519.2 


USE  OF  DECORRELATIVE  PROPERTIES  OF  DISCRETE  SPECTRAL  TRANSFORMATIONS  IN 
MULTI-ALTERNATIVE  RECOGNITION  OF  SIGNALS  IN  PRESENGE  OF  CORRELATED  NOISE 

Kiev  IZVESTIYA  VYSSHIKH  UCHEBNYKH  ZAVEDENIY:  RADIOELEKTRONIKA  in  Russian 
Vol  26  No  3,  Mar  83  (manuscript  received  after  revision  26  Apr  82)  pp  45-49 

[Article  by  A. I.  Rog  and  A. A.  Sirota] 

[Text]  A  method  is  presented  for  estimating  the  error 
probability  in  signal  recognition  against  the  background 
of  correlated  noise  with  incomplete  description  of  the 
correlation  connections  of  the  features.  The  proposed 
method  is  used  to  substantiate  the  possibility  of 
replacing  optimal  decorrelation  transformation  by 
partially  decorrelative  discrete  spectral  transformations. 


The  problems  which  arise  in  a  number  of  practical  applications ,  such  as 
automatic  recognition  of  pulsed  signals  of  various  forms:  S^(t)=m^(t)+N(t) , 

i=l,L,  where  m^(t)  —  determinate  signal  component;  n(t)  —  normally 

distributed  noise  with  correlation  function  a^R(T) ,  involve  using  decision 
rules  of  the  form  [1,2] 


^=?-inin[(;c  • 


•  nti)^ M  '  (x  ■ 


(1) 


In  (1)  x=(S(tQ) , . . . ,  S(t^_j^))  —  vector  of  amplitude  (feature)  samples  of 
signal  S(t)  obtained  in  the  time  interval  [0,Tl;  m^  —  vector  of  mathematical 
expectancies  of  features;  M  —  covariation  matrix  of  features. 

The  use  of  optimal  decision  rules  with  structure  (1)  involves  definite 
difficulties  which  arise  in  estimating  the  matrix  M  and  which  also 
require  a  very  great  deal  of  time  to  implement  the  matrix  transformations  in 
(1)  during  the  decision  making  stage.  It  becomes  necessary  to  examine  a 
special  class  of  quasi-optimal  recognition  algorithms  based  on  a  certain 
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"disregard"  of  the  correlation  between  features,  which  in  the  present 
case  is  manifested  in  that  a  matrix  K  which  is  slightly  different  from 
the  actual  matrix  M  is  employed  in  (1). 


Let  us  now  examine  the  losses  which  occur  when  this  class  of  quasi-optimal 
algorithms  is  employed.  If  the  error  probability  of  pairwise  separation  of 
signals  for  optimal  algorithm  (1)  is  defined  by  the  familiar  formular  [2] 


Pl-  =  (I)  (a,.y),  atj  =  0,5  {nii  —  ttijfM  ' —  ntj), 


(b(x:)  =  fe  "dir 
y2n} 


we  find  this  probability  for  the  quasi-optimal  algorithm  from  the 
distribution  parameters  of  the  logarithm  of  the  likelihood  ratio  of  the 
two  signals  [2]  with  matrix  K  substituted  instead  of  matrix  M: 

In  [/  5)]  '  (mi  (m,-  +  (m^  -  rnj) . 


It  can  be  shown  that 

Ml  {[In  /  (x)]/0ft}  =  ±  0,5  (mi  —  mjf  IT'  (ttii  —  mj);  k=-i,  /; 
D  {[In  I  {x)]/0„}  =  (m;  —  mjfK~'MlC^  (m;  —  mj), 


and  the  probability  of  error  in  pairwise  separation  in  this  case  is 

(mi  —  mjfK~^  (mj  —  mj) 


(nil  —  mj)^K  ^MK  '  (nii  —  mj). 


0  k 

The  relationship  between  the  probabilities  P_  and  is  determined  by 

the  relationship  between  the  quantities  ex. .  and  3. .  which  can  be  represented 
in  the  form  below  by  introducing  the  difference  vector  a=m^-m  ,  the  vector 


b=K  ^^^a  and  the  matrix  B=K 


aij  =  Q,^V'a^K  ^  a  =  o,5  'K, 

_±  _1_  fZ  _±  _  ‘  ZZirZZTZ 

^ij=Z*aZK  "K  ^a/2\/  a^K  'M/C  ' /C  ^  a=  b^bl2V  b^Bb . 
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In  order  to  estimate  the  losses  occurring  with  the  use  of  this  class  of 
quasi-optimal  algorithms  we  examine  the  quantity  A=a.  .£ot.  .  [l-min(3 .  . /ot  . )  ] 

for  a  given  error  probability  for  optimal  algorithm  and  corresponding 

value  of  a... 

13 

Using  Bekkenbakh  and  Kassel's  inequality  [3],  we  obtain 


1  <  a2./p2.  =  {b'^Bb)  {b^B  'b)i{b  bf  < - 


(Iq  + 


(2) 


where  —  minimum  and  maximum  eigenvalues  of  the  matrix 

g=|r-l/2j^-l/2^ 

We  obtain  from  (2)  the  following  formula  for  the  losses 
occurring  with  incomplete  allowance  for  feature  correlation  in  the 
quasi-optimal  recognition  algorithms: 

~  {<^u) ,  max  P'lj  =  O  {201,1  j  ( 3 ) 

V  ^0  +  ^A'-l  /  ' 


Thus,  for  an  algorithm  which  does  not  allow  fully  for  sample  correlation 
(K=I,  where  I  —  identity  matrix) 


ft  min  [(x  —  {x  —  m-^\  x^Oi 


71' 


(4) 


we  can  write  the  following  expression  for  the  maximum  possible  recognition 
error  probability  (p9,=const): 


P  7  =  a)  a,. 


2  y  X|,xv_i  \ 

><0  +  ^N-\ )  ' 


1 


(5) 


where  K^,  —  eigenvalues  of  matrix  M. 

The  difficulties  involved  in  using  information  about  the  correlation  of  the 
initial  features  —  signal  amplitude  samples  —  can  be  avoided  by  using 
the  familiar  decorrelation  method.  An  optimal  decorrelative  transformation 
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is  the  expansion  with  respect  to  discrete  eigenfunctions  of  the  matrix  M 
using  a  Karunen-Loew  expansion  [2] .  The  optimal  decision  rule^in  the 
space  of  the  new  uncorrelated  features  assigned  by  the  vector  y  appears  as 


k  min  [{y  —  ^  {ij  —  n^)] 


(6) 


where  y=Vx,  n^=Vm^,  V  —  optimal  decorrelative  transformation 

(ODT)  matrix.  The  main  shortcoming  of  the  ODT  is  the  difficulty  of  finding 
the  expansion  functions  and  the  excessive  amount  of  time  required  to 
execute  the  ODT.  It  is  therefore  better  to  try  to  replace  the  ODT  with 
some  other  orthogonal  transformation,  which  may  have  slightly  poorer 
decorrelative  properties  but  which  allows  the  transformation  procedure  to 
be  speeded  up  and  simplified  significantly.  Such  transformations  include 
Fourier,  Walsh  and  Haar  discrete  spectral  transformations  (DST) ,  for  which 
’*fast**  transformation  methods  have  been  worked  out.  The  idea  of  the  ODT -DST 
replacement  is  based  on  certain  facts  of  the  similarity  of  the  structure  of 
the  bases  of  orthogonal  functions  of  the  optimal  decorrelative  transformation 
and  of  spectral  transformations  assigned  on  a  discrete  finite  interval  [4,5]. 
For  for  example,  [5]  notes  that  the  Fourier  transform  approaches  the 

Karunen-Loew  transform  asymptotically. 

The  problem  of  estimating  the  effectiveness  of  the  use  of  an  approximate 
decorrelative  transformation  instead  of  the  ODT  is  equivalent  to  the  problem 
of  estimating  the  losses  which  occur  with  an  incomplete  description  of  the 
correlation  connections  of  the  signal  samples  in  quasi-optimal  recognition 
algorithms.  Actually,  replacing  the  ODT  with  some  DST  assigned  by  the 
matrix  H  actually  consists  of  using  in  place  of  optimal  algorithm  (6)  an 
analogous  structure  of  a  quasi-optimal  algoj^ithm  implemented  in  the  space 
of  the  spectral  coefficients  of  the  DST  —  z: 


k  min [(2  ^  {z-li)]^z€  Oft. 


(7) 


where  z=Hx,  E=diagP,  P=HMH^  --  covariation  matrix  of  spectral 

coefficients. 

Algorithm  (7)  does  not  make  complete  allowance  for  the  correlation  between 
the  spectral  coefficients,  and  is  quasi-optimal  with  respect  to  an  algorithm 
with  structure  (1);  however,  it  is  realized  in  a  spectral  coefficient 
space  in  which  the  covariation  matrix  of  these  coefficients  P  is  replaced 
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by  diagonal  dispersion  matrix  E.  We  can  therefore  estimate  the  losses  which 
occur  when  quasi-optimal  algorithm  (7)  is  used  instead  of  optimal  algorithms 
(1)  and  (6)  on  the  basis  of  expression  (3) ,  which  is  fairly  general.  The 
expression  for  the  maximum  possible  recognition  error  probability  for 
algorithm  (7)  for  a  given  value  of  —  the  error  probability  of  algorithm 

(6)  —  which  determines  the  effectiveness  of  the  ODT-DST  substitution,  looks 
like 


(““Sir)’ 


(8) 


-1/2  -1/2 

where  Xq»  —  eigenvalues  of  matrix  E  PE  '  . 


Fig.  1 


Formula  (8)  was  used  as  the  basis  for  investigating  the  effectiveness  of 
using  Fourier  and  Walsh  DST  in  quasi-op timal  signal  recognition  algorithms. 
Figure  1  shows  the  probabilities  of  pairwise  signal  recognition  errors  for 
different  quasi-optimal  algorithms  as  a  function  of  the  ratio  of  the  noise 
correlation  time  to  the  signal  duration.  For  comparison,  the  error  probability 
was  calculated  simultaneously  for  quasi-optimal  algorithm  (4) ,  which  does 

not  make  complete  allowance  for  sample  correlation.  Figure  1  a,  b  and  c  show 


P^  and  P. .  as  a  function  of  the  ratio  of  the  noise  correlation  time  to  the 
ij  ij 
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signal  duration  —  T,/T  for  different  values  of  N  and  types  of 

^  0 

correlation  function  R(t)  for  a  value  of  P  =0.01  (the  dotted  line  in 

^  J  — (^  'I 

Fig.  la  corresponds  to  N=32,  and  the  solid  line  to  N=l^  for  R(T)=e  ; 

In  Fig.  1  b  and  c  for  R(T)=e-aKl  cos  WT  and  R(T)=e"®'’^  ,  respectively). 

It  follows  from  these  relationships  (1 —  probability  P^|;  2F,  2W  — 

probability  computed  for  Fourier  and  Walsh  transformations,  respectively), 

that  the  losses  which  occur  with  the  use  of  quasi-optimal  algorithms  with 
structure  (7)  in  Fourier  and  Walsh  coefficient  space  are  significantly 
smaller  than  for  quasi-optimal  algorithms  (4)  realized  in  the  sample 
space.  The  decorrelative  properties  of  the  DST  are  particularly  strong 
for  small  values  of  T  /T,  when  the  quantity  P^^  is  practically  the  same 

as  P:..  As  follows,  e.g.,  from  Fig.  la,  by  comparison  with  algorithm  (6), 

based  on  the  use  of  the  ODT,  signal  recognition  accuracy  for  noise 
described  by  an  exponential  correlation  function  is  no  more  than  3-5% 
lower  when  quasi-optimal  algorithms  (7)  with  Fourier  transformation  are 
employed,  and  no  more  than  5-10%  lower  for  Walsh  transformation  for  values 
of  Tj^/T_<0.4. 
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APPROXIMATION  OF  FOURIER  TRANSFORM  KERNEL  WITH  HELP  OF  TWO-LEVEL  FUNCTIONS 
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Vol  26  No  3,  Mar  83  (manuscript  received  after  revision  5  Apr  82)  pp  57-59 

[Article  by  V.B.  Dmitriyev-Zdorov  and  M.N.  Surkov] 

[Text]  Spectrum  analysis  devices  can  be  simplified  significantly  by 
approximating  the  transformation  kernel  with  piecewise-constant  functions 
with  a  small  number  of  levels  [1].  It  is  particularly  attractive  to  do  the 
approximation  with  the  help  of  a  periodic  two-level  signal:  this  makes  it 
possible  to  increase  dynamic  range  significantly  and  to  simplify  the 
analyzer  structure.  The  spectrum  of  such  a  periodic  two-level  signal  must 
contain  a  first  harmonic,  after  which  the  series  of  harmonics  (2,3,...,p) 
is  lacking,  and  the  amplitudes  of  the  higher  harmonics  (p+1,  p+2,...,) 
are  arbitrary. 

The  objective  of  the  present  study  is  to  synthesize  a  periodic  digital 
signal  whose  spectrum  does  not  contain  harmonics  2,  3,  ...,  p.  We  shall 
represent  the  two-level  signal  which  has  the  required  spectral  structure 
in  the  form 


i=l 


,  ,  ,  .  rl»'r6  (0,t-i) 

where  h(T)={  ^ 

0,  for  other  x  —  describes  the  shape  of  rectangular  pulses 
with  duration  t^,  the  spacing  of  which  in  the  periodic  sequence  of 

pulses  y(t)  with  period  T=Nt^  is  determined  by  the  values  of  the  coefficients 


I,  =  {0.1}. 


(1) 


« 


29: 


We  represent  the  complex  amplitude  of  the  ith  harmonic  as 


T 

6 


/=1 


For  i<N  the  condition  is  equivalent  to  the  requirement 


N 

l^\ 


(2) 


The  simplest  signals  which  correspond  to  the  solutions  of  (2)  are  groups 
of  m  pulses  spaced  apart  by  N/(mi)  pulse  positions,  i.e.,  by  1/m  part  of 
the  period  of  the  ith  harmonic  component  (Fig,  1),  The  requirement  of 
(2)  is  satisfied  for  the  signal  shown  in  Fig.  1,  since  the  left  part  of 
(2)  can  be  viewed  as  the  sum  of  m  vectors  on  a  complex  plane  (Fig,  2)  with 
norm  of  unity  and  arguments  2Tm/m,  n=l,2,..,,  m  (which  can  be  obtained  by 
assuming  £=nN/ (mi)  in  (2)). 


Fig.  1 


Fig-  2 


We  note  that  the  condition  C^=0  is  equivalent  to  having  coefficients  a^  and 

b^  of  zero  in  the  expansion  of  y(t)  into  a  real  Fourier  series,  since  2C^= 

=a.-jb .  [2] . 

.1  1 
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The  other  solutions  of  (2)  can  be  obtained  by  linear  combination  of  the 
simple  solutions;  however,  (1)  must  not  be  violated  when  this  is  done. 


Let  be  a  group  of  m  pulses  corresponding  to  the  solutions  of  equation  (2) 

for  a  given  i.  In  a  periodic  signal  this  group  repeats  every  N  positions. 

In  the  spectrum  S^(2TTfQk)  of  such  a  periodic  signal  U^(t)  the  amplitude  of 

the  ith  harmonic  is  zero:  S^(27rfQi)=0,  where  fQ=l/T  —  frequency  of  first 

harmonic  of  signal  y(t).  The  convolution  spectrum  W(t)  of  such  signals  U^, 
i— 2 , . . .  ,p 


^  (0  =  ^2  (0  ^  (0  ^  ^  (0 


(3) 


p 

is  equal  to  the  product  of  the  spectra  of  these  signals  s  (2KfQk)  =  Y\  S.(2KfQk) 

1=2 

and,  consequently,  the  amplitudes  of  the  harmonic  components  i=2,...,p  in 
it  are  zero.  The  convolution  operation  in  (3)  is  the  discrete  convolution 

N 

of  periodic  pulse  sequences  of  the  type:  j=zi,...,N. 

i=\ 


The  signal  W(t)  is  also  periodic,  with  period  T=Nt^.  Obviously,  the  number 

of  positions  N  within  a  period  of  this  signal  is  a  multiple  of  all  factors 
of  the  type  m^A,  i=2,...,p. 

In  general  terms,  this  is  the  algorithm  for  constructing  the  sought  signal 
W(t).  It  should  be  noted  that  the  presence  in  the  convolution  (3)  of  a 
certain  group  U.,  as  will  be  shown  later,  makes  the  amplitude  of  the 
harmonic  numbered i  equal  to  zero,  as  well  as  a  series  of  higher  harmonics. 
This  makes  it  possible  to  reduce  the  number  of  groups  in  the  right  part  of 
the  convolution  (3)  and,  consequently,  to  reduce  significantly  the  number  of 
positions  N  in  the  period  of  the  signal  W(t),  which  is  extremely  important 
for  practical  purposes. 

Let  us  now  examine  the  spectral  structure  of  a  T-periodic  sequence  of  groups 
of  m  pulses  spaced  1^=1/ (im)  apart.  Finding  the  frequencies  of  the  harmonics 

which  are  lacking  in  the  spectrum  of  this  sequence  involves  solving  this 
equation: 


F  (f)  + ...  +  X  +  1)  =  0, 


where  X=e 


j  27rfT/  (j^  ) 
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and  corresponds  to  a  delay  operator  by  T/(iin),  while  F(f)  is  determined 
by  the  shape  of  the  pulses. 

Obviously,  the  equation 


+  X”'~^  +  ...  +  X  +  1  =(X"‘  —  I)/(X  —  1)  =  0. 


(4) 


is  a  special  case  of  (2)  and  has  the  solutions  k=l,2 , . . .  ,m-l. 

There  corresponds  to  each  of  the  solutiqns  of  (4)  a  series  of  frequencies 
which  are  found  from  the  formula  Xj^=x=ed^^^^' .  However, 

consequently,  2TTfT/ (im)=2TTk/m+27rr ,  r=0,  +1,  +2,...,  whence 

f  =  {k-\- rm)ilT,  A  =  1, 2 ,  m  —  l,  /•  _  o,  ^  2 


We  note  that  the  factor  (k+rm)  in  the  right  part  of  (5)  can  take  on  any 
integer  values  except  for  integer  multiples  of  m.  Consequently,  the 
presence  of  a  group  consisting  of  m  pulses  in  the  convolution  (3)  leads 

to  compensation  in  the  resulting  signal  W(t)  of  all  of  the  harmonics  whose 
numbers  are  multiples  of  1,  except  for  those  whose  frequencies  are 
represented  in  the  form  f^=fQlm£,  i!/=0,l,2 .  In  order  to  construct  a 

signal  W(t)  in  which  the  amplitudes  2,3,...,p  of  the  harmonics  are  zero 
it  is  necessary  to  perform  a  convolution  of  the  type  (3)  which  includes 
groups  of  m  pulses  such  that  i  runs  through  all  of  the  primes  and  all 
of  the  integer  powers  of  m  not  exceeding  p.  Let  I(p,m)  denote  the  set 
consisting  of  these  values  of  i.  As  was  noted  earlier,  the  number  of 
positions  N  within  a  period  of  the  signal  W(t)  must  be  a  multiple  of  all 
factors  of  the  type  im.  Therefore,  the  minimum  N  is  defined  as  N=mLCM 
{l(p,m)},  where  LCM  —  least  common  multiple. 

As  m  increases  the  number  of  groups  in  convolution  (3)  becomes  smaller. 

However,  the  resultant  signal  W(t)  remains  a  two-level  signal  only  for 
m=2  and  m=3 . 

The  highest  suppressed  frequency  fj^=fQp  is  in  principle  unbounded.  However, 

for  large  p(p>10)  the  required  number  of  positions  N  in  a  signal  period 
increases  immeasurably. 

Table  1  shows  the  minimum  N=N(p)  for  the  case  in  which  the  signal  W  (t) 

^13 

is  synthesized  by  convoluting  simple  two-  and  three-pulse  signals  ^ 
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(for  m=TQ.^=2  and  iii=in2=3)  and  q  —  number  of  pulses  in  period  of  W^(t). 


TABLE  1 


p 

\nu/ni^ 

1  ‘ 

NJN, 

P 

NJN^ 

2 

2/3 

4/6 

2/3 

9 

-/3 

1680/1890 

64/243 

3 

2/3 

12/18 

4/9 

10 

— / — 

1680/1890 

64/243 

4 

2/- 

24/18 

8/9 

11 

2/3 

18480/20790 

128/729 

5 

2/3 

120/90 

16/27 

12 

— / — 

18480/20790 

128/729 

G 

— / — 

120/90 

16/27 

13 

2/3 

240240/270270 

256/2187 

7 

2/3 

840/630 

32/81 

14 

— 

240240/270270 

256/2187 

8 

2/- 

1680/630 

64/81 

15 

— / — 

240240/270270 

256/2187 

16 

2/- 

480480/270270 

512/2187 

The  class  of  two-level  signals  W^(t)  can 


be  used  for  spectrum  analysis 


of  signals  occupying  a  limited  frequency  band  f  ,  <f<f  .  It  can  be 

shown  that  there  is  no  methodical  analytical  error  if  f^^.^-1, 

which  confirms  the  advisability  of  adapting  the  synthesized  signal 

to  the  analyzed  process.  One  feature  of  this  approach,  in  contrast  to 
universal  algorithms  (specifically  the  FFT  [3,4]  where  the  transformation 
kernel  is  not  approximated)  is  the  fact  that  by  increasing  the  number  of 
discrete  samples  of  the  analyzed  signal  (increasing  the  digitization 
frequency)  it  is  possible  to  eliminate  the  multiplication  operation 
completely.  In  order  to  obtain  the  expansion  coefficients  of  the 
analyzed  signal  it  is  sufficient  to  add  the  values  of  the  samples  taken 
at  the  instance  corresponding  to  non-zero  levels  of  the  digital  signal 
W  (t).  In  calculating  different  expansion  coefficients  of  the  signal 
P 

in  question,  the  samples  can  be  added  Independently  and  in  parallel. 


Increased  speeds  can  be  expected  in  analyzing  signals  with  relatively 
narrow  bandwidths,  especially  when  the  code  word  lengths  in  the  computer 
are  long.  It  should  also  be  pointed  out  that  it  is  possible  in  principle 
to  filter  the  analyzed  signal  first  in  order  to  represent  it  as  the  sum 
of  narrowband  components,  and  then  to  apply  the  analysis  algorithm  in  question. 
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EFFECTIVENESS  OF  DIGITAL  TRANSMISSION  EMPLOYING  MODIFIED  WAGNER  METHOD 
TO  CORRECT  SINGLE  ERRORS 
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Vol  26  No  3,  Mar  83  (manuscript  received  after  revision  23  Sep  82)  pp  70-72 

[Article  by  R.E.  Gut,  M.Ya.  Lesman  and  A.S.  Shutov] 

[Text]  Various  decoding  procedures  which  allow  for  the  reliability  of  the 
decisions  made  about  the  signal  elements  transmitted  are  employed 
extensively  in  digital  transmission  systems.  For  transmission  over  radio 
channels  these  procedures  are  especially  convenient  when  multi-channel 
modems  are  used  (such  as  the  Kineplex,  MS  [expansion  not  given],  etc.  [1]). 
In  contrast  to  the  approach  usually  taken  [2],  it  is  best  to  judge  the 
symbol  reliability  according  to  the  estimates  of  the  signal/noise  ratio 
provided  by  meters  which  are  specially  added  to  each  of  the  modem  channels. 

The  present  study  solves  the  problem  of  estimating  the  noise  tolerance  of 
parallel  reception  of  signals  coded  with  an  (M,  M-I)  code  and  decoded  by 
the  Wagner  method.  The  indicators  which  determine  noise  tolerance  are  the 
probabilities  of  correct  reception  of  a  code  combination  and  a  symbol  in  an 
arbitrary  phased  channel.  Channel  fading  is  assumed  to  be  statistically 
uniform  and  independent. 

It  can  be  demonstrated  that  the  expressions  for  the  mean  probahilltles  of 
correct  reception  of  a  code  combination  and  of  a  symbol  in  an  arbitrary 

jth  channel  in  the  case  in  question  appear  as 


Qi  =  (l 


_  “ "  f” 

+  M  J  J  p(h)f  (S,  /i)  d/i|  r  (i_p(/,)]  F  (S\h)dh 
0  0  lo-' 


M—l 


dS. 


(1) 


where  p  —  mean  error  probability  during' reception  of  elementary  signal 
in  any  of  M  statistical  uniform  channels ;  p (h)  —  current  error  probability 
in  receiving  elementary  signal  as  function  of  actual  signal/noise  ratio; 

M  — —  number  of  symbols  in  code  word  (number  of  parallel  channels) ;  f (h) 
density  function  of  signal/noise  ratios;  f(S,h) —  combined  density 
function  of  estimate  of  signal/noise  ratio  and  actual  value  of  slgnal/nolse 
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ratio;  F(S|h)  —  conditional  distribution  function  of  signal/noise  ratio 
estimate  assuming  that  the  actual  signal/noise  ratio  is  h^; 


00  W  ^ J 

<32  =  (1-p)-4-  f  /(5)|^'^“‘(S)-[f(S)-2  {p(h)f(h)F(S\h)dh\  [dS  + 

0  ^0  (2) 
CO  00 

+  i[  (S)  ^  p  (h)  f  (S,  h)  dhdS, 

6  6 


where  f(S)  and  F(S)  —  unconditional  density  function  and  distribution 
function  of  estimate  of  current  value  of  signal/noise  ratio,  respectively. 


0  ^  8  12  16  n 


) 

. 

<5- 

Fig.  1  Fig.  2 

Let  us  now  specify  the  formulas  derived  above  for  the  case  of  optimal 
incoherent  reception  of  signals  employing  single  phase-differential  keying 
[2].  We  shall  assume  that  the  channel  fading  is  described  by  a  Rayleigh 
distribution,  and  that  the  receiving  location  employs  the  devices  for 
measuring  the  current  signal/noise  ratio  which  are  described  in  [3],  for 
which 


/  (S  I  ft)  =  [1/r  (v)]  (4vft=)'’S'^-'  exp  (—  4vft2S). 


(3) 
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where  2v  —  size  of  sample  used  to  form  estimate  of  signal/noise  ratio. 

By  substituting  (3)  in  (1)  and  (2),  we  obtain  the  following  analytical 
formulas : 


Qi 


1  '+1  r+^i 

^  x'^^dx  * 

[i_.  .  1  ] 

l2(l+/i2)J  2., 

(4) 


M-l 


+  2/1^ 


2(1+/i2)  2M  +  2  '[*  — 

+  -5-1 


nM-1 

^vM~\dx 


(5) 


Calculations  were  done  in  accordance  with  (4)  and  (5),  some  of  the 
results  of  which  are  shown  in  Figs.  1  and  2.  Figure  1  shows  (1-Q^) 

during  reception  of  a  code  combination  as  a  function  of  the  sample  size 

(n“2v)  for  M=2  (dotted  lines)  and  M=6  (solid  line)  and  for  two  values  of 

h^  =10  (curves  1  and  3)  and  h„^=100  (curves  2  and  4) .  Figure  2  shows 

^  ^2  2 
the  same  for  (I-Q2)  with  M=6  and  h^  =10  (curve  5)  and  h^  =100  (curve  6). 

We  note  that  when  m=2  the  probabilities  and  are  the  same.  Furthermore, 

the  straight  dot-and-dash  lines  in  Figures  1  and  2  show  the  values  of  the 
corresponding  error  probabilities  during  code  combination  and  symbol 
reception  which  are  achieved  in  a  hypothetical  system  belonging  to  the  class 
in  question  in  which  precise  information  is  available  about  the  actual 
signal/noise  ratios  in  each  channel. 

Analysis  of  the  findings  permits  the  following  conclusions: 

1.  The  use  of  the  modified  Wagner  procedure  (considering  the  signal/noise 
ratio  meters)  allows  the  error  probability  for  code  combination  and  symbol 
reception  in  an  arbitrary  channel  to  be  reduced  by  almost  an  order  of 
magnitude,  even  in  statistically  uniform  channels.  This  procedure  should 
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provide  even  greater  effectiveness  when  one  of  the  channels  is  affected 
by  spectrally  concentrated  noise. 

2.  The  mean  error  probabilities  for  code  combination  (1-Q^)  and  symbol 
(l-Q^)  reception  in  an  arbitrary  channel  are  practically  the  same  as  the 

corresponding  figures  for  the  hypothetical  system  with  a  sample  size  of 

n=2  =15-20. 
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